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__3 \ Abstract 

^p H " We present a phenomenological study of the current status of the application of our 
Oh- approach of exact amplitude-based resummation in quantum field theory to precision 
QCD calculations, by realistic MC event generator methods, as needed for precision LHC 
l— ~~ '! physics. We discuss recent results as they relate to the interplay of the attendant IR- 
Improved DGLAP-CS theory of one of us and the precision of exact NLO matrix-element 
matched parton shower MC's in the Herwig6.5 environment as determined by comparison 
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^Nj ■ to recent LHC experimental observations on single heavy gauge boson production and 
decay. The level of agreement between the new theory and the data continues to be a 
reason for optimism. In the spirit of completeness, we discuss as well other approaches 
to the same theoretical predictions that we make here from the standpoint of physical 
precision with an eye toward the (sub-)l% QCDcg)EW total theoretical precision regime 
for LHC physics. 
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1 Introduction 



With the recent announcement [T] of an Brout-Englert-Higgs (BEH) |2j candidate boson 
after the start-up and successful running of the LHC for three years, we have entered 
the era of precision QCD, by which we mean predictions for QCD processes at the total 
precision tag of 1% or better. Given the expected role of precision comparison between 
theory and experiment in determining the detailed properties of the newly discovered 
EBH boson candidate, the attendant need for exact, amplitude-based resummation of 
large higher order effects is now more paramount. We have argued in Refs. [31H] that such 
resummation allows one to have better than 1% theoretical precision as a realistic goal in 
such comparisons, so that one can indeed distinguish new physics(NP) from higher order 
SM processes and can distinguish different models of new physics from one another as 
well. In what follows, with an eye toward its interplay with NLO exact, matrix element 
matched parton shower MC precision issues, we present the status of this approach to 
precision QCD for the LHC in connection with its attendant IR-improved DGLAP-CS [5] 
E] theory 018] realization via HERWIRI1.031 [9] in the HERWIG6.5 [Hi] environment. 
We realize here the attendant exact, NLO matrix element matched parton shower MC 
realizations for both HERWIRI1.031 and HERWIG6.5 in the MC@NLO p] methodology, 
wherewith we make the corresponding comparisons with recent LHC data that we present 
herein. 

In this way, we continue the strategy of building on existing platforms to develop and 
realize a path toward precision QCD for the physics of the LHC. What we exhibit is an 
explicit union of the new IR-improved DGLAP-CS theory and the MC@NLO realization 
of exact NLO matrix element (ME) matched parton shower MC theory. We note that 
we are also pursuing the implementation p2] of the new IR-improved DGLAP-CS theory 
for HERWIG++ [13], HERWIRI++, for PYTHIA8 [H] and for SHERPA [T5], as well as 
the corresponding NLO ME/parton shower matching realizations in the POWHEG [16] 
framework. For, our ultimate goal is a provable precision tag on our theoretical predictions 
and one of the strongest cross checks on theoretical precision is the difference between two 
independent realizations of the attendant theoretical calculation. We will present such 
cross checks elsewhere [T2] , 

The paper is organized as follows. We motivate the development of the new IR- 
improved DGLAP-CS theory in our approach to precision QCD theory by showing in 
the next section how it arises naturally in the effort to realize a provable precision from 
our approach [I] to precision LHC physics. This will allow us to expose properly the 
interplay between the NLO ME matched parton shower MC precision and the new IR- 
improved DGLAP-CS theory. In the interest of completeness, we also review in the next 
section the relevant aspects of the approach in Ref. [3], which is an amplitude-based 
QEDg)QCD(= QCD ® QED) exact resummation theory realized by MC methods. This 
paper is an extension of a shorter report of our results and analysis which is given in 
Ref. [17J. Accordingly, also in the interest of completeness, we discuss as well in the next 
section other possible approaches to the predictions we make herein from the standpoint 
of their physical precision. We believe this discussion will aid the reader in putting the 
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results we present in a proper context. Section 3 contains the applications to the recent 
LHC data on single heavy gauge boson production with an eye toward the analysis in 
Refs. [9] of the analogous processes at the Tevatron. In this paper we will focus on 
the single Z/7* production and decay to lepton pairs for definiteness. The other heavy 
gauge boson processes will be taken up elsewhere [12] . Section 4 contains our concluding 
observations and our outlook. 

2 Review of Our Approach to Precision LHC Physics 

We take the starting point for what we discuss here as the following fully differential 
representation of a hard LHC scattering process: 

da — / dx 1 dx 2 F i (x 1 )F j (x 2 )da res (xiX2s), (1) 

where the {Fj} and da ves are the respective parton densities and reduced hard differential 
cross section and where we use the subscript on the latter to indicate that it has been 
resummed for all large EW and QCD higher order corrections in a manner consistent with 
achieving a total precision tag of 1% or better for the total theoretical precision of ([1]). The 
key theoretical issue for precision QCD for the LHC is then The proof of the correctness 
of the value of the total theoretical precision Aa t h of ([I]) is then the key theoretical issue 
for the realization of precision QCD for the LHC. The theoretical precision in question 
can be represented as follows: 

Aa th = AF © Aa ves (2) 

where AA is the contribution of the uncertainty on the quantity A to Aoth- Here, we 
stress that we discuss the situation in which the two errors in (T5]) are independent for 
definiteness; (T5]) has to be modified accordingly when they are not. We further stress 
that the proof of the correctness of the value of the total theoretical precision Aa t ^ 
is indispensible in order to validate the application of a given theoretical prediction to 
precision experimental observations, for the discussion of the signals and the backgrounds 
for both Standard Model(SM) and new physics (NP) studies, and more specifically for 
the overall normalization of the cross sections in such studies. If one uses a calculation 
with an unknown value of A<r t h for the attendant studies, the NP itself can be missed. 
We do feel a need to continue to state here that this point simply cannot be emphasized 
too much. 

In the interest of completeness, we note that, by our definition, which follows the 
discussion in Ref. [TS], Acr t h is the total theoretical uncertainty that derives from the 
physical and technical precision contributions: the physical precision contribution, Aaf^ yB , 
arises from such sources as missing graphs, approximations to graphs, truncations,....; 
the technical precision contribution, Ao"{£ ch , arises from such sources as bugs in codes, 
numerical rounding errors, convergence issues, etc. We want to observe that here, when 
we reference bugs in codes we have in mind that all gross errors such as those that give 
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obviously wrong results, as determined by cross checks, are eliminated and we have left 
programming errors such as those in the logic: suppose for programming error reasons 
a DO-loop ends at 999 steps instead of the intended 1000 steps, resulting in a per mille 
level error, that could alternate in sign from event to event. As per mille level accuracy is 
good enough in many applications, the program would remain reliable, but it would have 
what we call a technical precision error at the per mille level. With this understanding, 
the total theoretical error is then given by 

Aa th = A< hys © A<<\ (3) 

Although the desired value for Aoth depends on the specific requirements of the observa- 
tions, as a general rule, it should fulfill Acr t h < /Acr expt , where Acr exp t is the respective 
experimental error and / < ^, so that the theoretical uncertainty does not significantly 
adversely affect the attendant physics analysis of the data. 

In the interest of completeness, we will now recapitulate the theory we have developed 
to realize such precision rigorously. This theory is the QCD <g> QED resummation theory 
in Refs. [!] for the reduced cross section in (TjQ) and for the resummation of the evolution 
of the parton densities therein as well. More precisely, the master formula, which applies 
both to the resummation of the reduced cross section and to that of the evolution of the 
parton densities, may be identified as 

_ „SUM m (QCED) y-oo _j_ r ™ d 3 k H 
uu res — e Z-m,m=0 n \m\ J llj'i=l k h 

11?2=1 k' j2 J (27r) 4Ci 

Pn,m{kl, ■ ■ ■ , k n ] k[, . . . , ^m)~^jp^p> (^) 

where da ves is either the reduced cross section da res or the differential rate associated 
to a DGLAP-CS [51E] kernel involved in the evolution of the {Fj} and where the new 

(YFS-style [13120]) non-Abelian residuals j3 ntTn (ki, k[, ... , k' m ) have n hard gluons 

and m hard photons and we show the final state with two hard final partons with mo- 
menta P2, q_2 specified for a generic 2f final state for definiteness. The infrared functions 
SUM IR (QCED), Dqced are defined in Refs. |HE1IB] as follows: 

SUM IR (QCED) = 2aMB^ ED + 2a a B$ ED 

/dk 
- (e-^ - 9(K max - k )) S$ ED (5) 

where the dummy parameter K max is such that nothing depends on it and where we have 
introduced 

jjnls Tjnls I rynls 

&QCED — D QCD ~r &QEDi 
Tynls _ Tjnls , r>nls 
QCED — n QCD ~r D QEDi 

nnls qtlIs , onls fr?\ 

D QCED — '-'QCD i '-'QED- \°) 

Here, the superscript nls denotes that the infrared functions are DGLAP-CS synthesized 
as explained in Refs. [4f7|l8|l2l] and the infrared functions Ba, Ba, Sa, A = QCD, QED, 
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are given in Refs. JUCOEJESED]- We stress that simultaneous resummation of QED and 
QCD large IR effects is exact here. 

We may describe the physical meanings of the key components in the master formula 
(BJ as follows. In the language of Ref. [22J, the exponent SUMtr(QCED) sums up to the 
infinite order the maximal leading IR singular terms in the cross section for soft emission 
below a dummy parameter K majX and the exponent -Dqced does the same for the regime 
above K max so that (jlj) is independent of K max - it cancels between SUM IR (QCED) and 
-DqcedIE I n order to maintain exactness order by order in perturbation theory in both a 
and a s in the presence of these resummed terms, as explained in Refs. [U0IE] we generate 
the residuals fi n ,m by iterative computation to match the attendant exact results to all 
orders in a and a s . We need to emphasize the following most important point. In our 
formulation in (jlj) the entire soft gluon phase space is included in the representation - no 
part of it is dropped. 

As it is explained in Refs. [I], the new non-Abelian residuals f3 m>n allow rigorous 
shower/ME matching via their shower subtracted analogs: in (jlj) we make the replace- 
ments 

where the /3 njm have had all effects in the showers associated to the {Fj} removed from 

them. To see how we make contact between the fi n ,m and the differential distributions in 
MC@NLO we proceed as follows. We represent the MC@NLO differential cross section 
via [EE] 

do MC&NLO 



B + V + j{R MC -C)d® R d$ B [A MC {0) + J{R M c/B)A MC {k T )d$ R } 

+ (R- RMc)AMc(k T )d<S> B d<S>R 

where B is Born distribution, V is the regularized virtual contribution, C is the corre- 
sponding counter-term required at exact NLO, R is the respective exact real emission 
distribution for exact NLO, Rmc — Rmc{Pab) is the parton shower real emission distri- 
bution so that the Sudakov form factor is 

A MC (p T ) = e i-f^n RMc( i BA ' R) 0(kT^B^ R )- P T)] 

, where as usual it describes the respective no-emission probability. The respective Born 
and real emission differential phase spaces are denoted by d$A, A — B, R, respectively. 
We may note further that the representation of the differential distribution for MC@NLO 
in ([HD is an explicit realization of the compensation between real and virtual divergent 
soft effects discussed in the Appendices of Refs. [H[8] in establishing the validity of @ for 



3 If we want to include more of the maximal exponentiating terms from the formalism of Ref. 22 in 
the two exponents SUMir,(QCED), -Dqced, we may do so with a consequent change in the attendant 
residuals fl n ,m- 
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QCD - all of the terms on the RHS of (jHD are infrared finite. Indeed, from comparison 
with (TjJ restricted to its QCD aspect we get the identifications, accurate to 0(a s ), 



1 s 
2 

2^i,o — R — Raw — BSqcd 



o ,o = B + (B/A MC (0)) [ (RMc/B)A M c(k T )d$ R 

3 (8) 



where we defined ITT 



B = B(1- 2a s MB QCD ) + V + J (R MC - C)d<$> R 

and we understand here that the DGLAP-CS kernels in Rmc are to be taken as the 
IR- improved ones as we exhibit below 011]. Here we have written the QCD virtual 
and real infrared functions Bq C d and Sqcd respectively without the superscript nls for 
simplicity of notation and they are understood to be DGLAP-CS synthesized as explained 
in Refs. [UCOIE] so that we avoid doubling counting of effects. We also re-emphasize that 
we do not drop any effects here in We observe further that, in view of (jHJ), the way 
to the extension of frameworks such as MC@NLO to exact higher orders in {a s , a} is 

therefore open via our /3 n>m and will be taken up elsewhere [T2] . 

We point out that in Refs. [7H9] our methods for resummation of the QCD theory 
have been shown to be fully consistent with the methods in Refs. [23| [2l]. Specifically, it 
is shown in Refs. [7H9] that the methods in Refs. [23l|2l] give approximations to our hard 

gluon residuals /3 nj0 ; f° r > unlike the master formula in (JJJ), the methods in Refs. [23||2T] 
are not exact results. To see this, observe that the threshold-resummation methods in 
Refs. [23], using the result that, for any function f(z), 



i 

n-l 



dzz n - L f(z) 







<(-) max .|/(*)|, 
n ze[o,i] 



drop non-singular contributions to the cross section at z — > 1 in resumming the logs in 
n-Mellin space. The SCET theory in Refs. [21] drops terms of O(X) at the level of the 
amplitude, where A = yjK/Q for a process with the hard scale Q with A ~ .3GeV so that, 
for Q ~ lOOGeV, we have A = 5.5%. From the known equivalence of the two approaches, 
the errors in the threshold resummation must be similar. Evidently, we can only use 
these approaches as a guide to our new non-Abelian residuals as we develop results for 
the (sub-)l% precision regime. 

In view of the specific processes which we consider in the next Section, we will continue 
here with our consistency discussion as it relates to the theory of QCD resummation for 
the specific heavy gauge boson production-type processes offered in Refs. J25j[26] , which 
is again an approximate formalism that could be used to make approximations to our 

hard gluon residuals /3 nj0 . The theory in Refs. [25l[26] is in wide use at the LHC and 
in the data analyses for the Tevatron - see for example the recent analyses in Refs. 
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where this theory, as it is implemented in the MC integration program ResBos [2"8H5U] is 
compared to recent LHC data and to recent analyses of Tevatron data. We need to note 
that the theory in Refs. p5|26] builds on and extends beyond the considerable literature in 
Refs. (31] aimed at the analogous processes to those under study here. Could we perhaps 
employ this formalism in Refs. [2"5|l2"6] as well to reach < 1% physical precision QCD 
predictions? Let us note that the authors in Ref. [26J give us a hint to the answer to our 
question in their footnote on the journal page with number 215 for that paper, wherein 
they equate as 'negligible' 20% in discussing possible nonperturbative contributions in 
their formalism. Let us keep this footnote in mind in what follows. 

We first note that a defining formula for the approach in Refs. [251126] is that for the 
differential cross section for the pt distribution for the heavy gauge boson production 
in hadron-hadron collisions, which we specialize here to the case of the Drell-Yan 7* 
production for definitenes^], where we record the result in Eq.(2.2) from Ref. [26] : 

1CW T ~ ^{S m^ T ^ Q - XA - XB) + ! '«S^,«,«>} (9) 

where we have the usual kinematics so that Qt = Pt is the 7* transverse momentum, 
A,B are protons at the LHC, s is the cms squared energy of the protons, Q M is the 
7* 4-momentum so that Q 2 is the 7* mass squared, and y = |ln(Q + /Q") is the 7* 
rapidity so that xa = e y Q/y/s and xb = e~ y Q/y/s. We have in mind that Q is near 
Mz here. In ([9]), the term involving W carries the effects from QCD resummation as 
developed in Refs. [251126] and the Y term includes those contributions which are 'regular' 
at Qt = Pt — > in the sense of Refs. [25l|26], i.e., order by order in perturbation theory 
they are derived from the part of the attendant hard scattering coefficient that are less 
singular than Q^ 2 x (logs or 1) or 5(Qt) as Qt = Pt 0. Our question concerns the 
physical precision of the W term; for, the Y term is perturbative and can be computed 
in principle to thej^equired accuracy by the standard methods. 

The result for W given in Refs. [25|[26] is as follows. When we have b << 1/A where A 
is aj^ypical QCD hadronic mass scale such as the inverse of the proton radius, the result 
for W is [25l[26l 




x ^e^CjaixA/Uib; C 1 /C 2 ; g(fi), ^)C jb (x B /^B,b; C 1 /C 2 ; g(ji), fi), (10) 

j 

where we define [25] [26] 

S sc = ^lH^f)A(g(fi); Ci) + &, C 2 )} (11) 

4 The analogous results for the W ± and Z/j* are obtained by straightforward substitutional manip- 
ulations of the EW aspects of the formula we record here as described in Ref. [26] so that we omit such 
manipulations here without loss of content of the QCD aspects our discussion. 
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Here we show explicitly the dependences of the parton density functions {Fg} (note that 
our Fi(x) corresponds to fg(x)/x in Refs. [25l[26]), and we show as well the dependences 
of the perturbatively calculable exponentiation and scattering coefficient functions A, B, 
and C as they are given in Refs. [25]|26]. The scale /x is usually set to C\jb and C% and 
Ci are order 1 constants chosen to optimize the resultant perturbation expansions - see 
Refs. [25][26] for more discussion on this point. The authors in Ref [25]|26] note that this 
result for W(b), when 1/Q « b « 1/A, is accurate up to terms 0(m q /Q, \/{bQ)) 
where m q represents the quark masses. In addition, the latter authors argue that in the 
result OH]) they have dropped terms of 0(m q /Q, Qt/Q) in the regime where < Qt -C Q. 
These last two statements are seen explicitly in Sect. 9 of Ref. [25J in eq.(9.1) and the 
equation immediately preceeding it, for reference. 

When we have to consider as well the regime b > 1/A, the authors in Ref. [25] 126] 
argue that we have to replace the functions A, B, and C according to 

A{g{p)\ Ci) -»■ A{g{p),m q /pL] C x ) 
B(g(jj); C h C 2 ) -> B(g(jj),m q /jj; C u C 2 ) 

^2 d U I d£ B F a/A (U; v) F b/B{Z,B\ ^)^2e 2 j C ja (x A /U,b;C 1 /C 2 ;g(fx), n) 

a,b Jx A Jx b j 

x C^xs/iB^-d/Ci-giii)^) -> ^e J 2 P jM (x^,6;C 1 /C 2 )P j/B (a; B ,6;C , 1 /C 2 ) (12) 

j 

where we see that we can no longer neglect the quark masses and that we no longer 
have the convolutions C * F: the functions Vj/A reduce to the latter form when we take 
b « 1/A and drop terms 0(bA), according to the arguments in Refs. [2"5|l2l)]. 

The authors in Refs. [25||26] then join these two results, that for the regime b << 1/A 
and that for the regime b > 1/A, via an ansatz as follows. Defining 

W(b;Q,x A ,x B ) = Y^^sM^mxb), (13) 

3 

which we see from fflUj) and f fl2|) is well-defined for both regimes of interest in b, the 
authors in Ref. [26] use 

b* = Vv / TTWC^ (14) 

where 6 max is a parameter still in the perturbative regime where f flO]) holds, to write the 
joining ansatz 

Wj(b;Q,x A ,x B ) = W 3 (b*;Q,x A ,x B )e { - HQ2/Q ^ 9l{b) ~ 9 ^ A{xA ^ (15) 

where the functions g\, gj/ A , gj/ B have non-perturbative content and must be determined 
from data. With this ansatz, the authors in Ref. [26] arrive at the representation 

sew - ^{ / £^z&*™.™y*-^«^>^>^ m 



3 



+ Y(Q T ;Q,x A ,x B ) }. (16) 
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We point-out for completeness and illustration that in ResBos in Refs. [2"9~tl3"U] three 
realizations of the non-perturbative functions in f fl6|) were considered as follows: 



\n(Q 2 /Ql)g 1 (b) + g j/A (x A ,b) + g j/B (x B ,b) 

([ gi +g 2 HQ/(2Q ))]b 2 , DWS [32] 
= I [ 9l + g 2 \n(Q/(2Q ))}b 2 + gi g 3 \n(100x A x B )b, LY [29] 
{[g 1 +g 2 HQ/(2Q ))+g 1 g 3 \n(100x A x B )]b 2 , BLNY [30]. 



(17) 



We observe that the $ are parameters on the RHS of (117j) and that the best fit to the 
data considered in Ref. [30] was obtained from the BLNY parametrization with the values 




when 6 max = 0.5GeV~ and Qo = 1.6GeV. 

On the question of the physical precision of (ITS]) , we first observe that, with the 
errors shown on the constants gi for the best BLNY parmetrization, the —0.02 error on 
g 2 represents already a 1.5% uncertainty at the saddle point position for the integration 
over b in the respective Z production analogue of ([9]) found in Ref. [26] and recapitulated 
below in ( TT9l - this alone exceeds the theory error budget in precision QCD theory that 
we advocate here so that one would need these errors reduced considerably for them to 
be useful for our purposes. But, what is more important is that, on the LHS of (jTTj) 
the gi which do not (do) multiply ln(Q/Qo) are unspecified functions of Xj,b (b) that 
are required to vanish at b = while on the RHS these dependences are simplified to 
second order polynomials in b and either no or linear ln(x A x B ) dependence on xf these 
simplifications are generically ad hoc and can not be considered as a rigorous platform 
for testing the fundamental QCD theory. They(the simplifications) do not follow from 
the formalism of Ref. [251126] . As the authors in Ref. [32] have emphasized, the latter 
formalism could only be used to give a prediction for perturbative QCD in the regime 
where the details of the non-perturbative parametrization are not important. We can not 
emphasize this point too much. 

On the question of the physical precision of ( 1T61) we additionally note that in the regime 
< Qt "C Q it has an error of 0{Qt/Q) according to the authors in Ref. [25l[26], as 
suchjterms have been dropped, order by order in perturbation theory in the construction 
of Wj - we make reference again specifically to Eq.(9.1) and the equation immediately 
preceding it in Ref. [25J. For Q = Mz and Qt = 5GeV, this is an error of = 5.5% and its 
twice this size at Qt = lOGeV. This shows that the formalism in Ref. [2S1I2S] can not be 
used for (sub-)l% theory predictions for the heavy gauge boson pt spectrum at the LHC. 
One might think that as Qt —> 0, where these dropped terms should vanish, the theory 
in Refs. [25j[26] would be accurate enough for (sub-)l% precision QCD predictions. This 
is not true as well because at Qt = 0, following Parisi and Petronzio in Refs. [31J, the 
authors in Ref. [26] show that the integral over b in ([9]) is dominated by a saddle point at 




g x = 0.21^GeV 2 , g 2 = O.eS^GeV 2 , g 3 = -0.6 



+0.05 
'-0.04 



(18) 




(19) 



where we have A = J2n ^ N \ a s/^) N and f3\ = ^(33 — 2N F ) where N F is the number of 
effective quark ffavors, which we take as 5 here. Thus, we have [251126] A^ = 4/3 so that 



1 (Q- ~ - 41 



USP ~ A VA 

For Q = M z , we get 6 S p = .48GeV _1 when we follow Ref. [26] and take A = 0.15GeV. 
This result for bsp is in the perturbative regime. Moreover, we know from Refs. [25|l26] 
that the error on the Z production analogue of Wj at this point is 0(l/(bspMz)) = 2.3% 
- see again Eq.(9.1) in Ref. [25]. Thus, even when the perturbative regime obtains so 
that the result (TTU|) should be reliable, the error in it is too large for use as anything but 

a guide in constructing our precision QCD ® EW theory residuals /3 n>m in our master 
formula (j4]) for < 1% precision LHC physics. The results reported in Refs. [27] on the 
comparison of ResBos with recent ATLAS and Tevatron data related to the p F spectrum 
in single Z/j* are consistent with our estimates on the physical precision of ( lT6i) : ResBos 
misses the data by 2% for the p F near and it misses the data by > 5% for the regime 
of p T > lOGeV. 

The non-perturbative factor in (fT6j) compromises the predictive power of the formalism 
in Refs. [2SJ[2n]. The authors in Ref. [32] have argued that for W, Z production as 
considered here the prediction in f ll6p in the regime 6GeV < Qt < 16GeV is insensitive 
to the non-perturbative parametrization in f|T7|) . This leaves the effects we have discussed 
above as the main obstacles to using the theory in Refs. [25l,[26] for precision QCD physics 
predictions at the LHC in this latter regime when that precision is at or below 1%. 

We note that in Ref. [33], another version of ([9]) is presented which has the same 
physical precision limits as that realized in ResBos with a different treatment of the non- 
perturbative regime. As the comparison with the data shown in Refs. [27] shows in the 
perturbative regime, our estimates of this physical precision for this regime also apply to 
the realization of ([9]) in Ref. [33]: it cannot be used for 1% precision LHC physics studies. 

Here, let us note that the authors in Ref. [33J discuss the uncertainty in their results 
as a function of the variation of the (perturbative scale) values associated with their 
renormalization, factorization and resummation scales. For example, they estimate that 
this is at the level of 10% near the peak of the Qt spectrum in the single Z/j* production 
at the LHC. These uncertainties can in principle be reduced by computing the perturbative 
terms in their results to higher and higher orders. They also estimate an additional PDF 
error at this region at the level of 2%, which in principle can be reduced by improving the 
determination of the respective PDF's. We stress that the error in the defining result 
derived in Refs. [25j[26] on which their results are based that we discuss here is an error 
of 0(Qt/Q) that applies order by order in the perturbation theory - it is separate from 
the scale and PDF errors discussed in Ref. [33] . 

Indeed, in Ref. [31], the SCET approach is used to recover Q and it is shown explicitly 
that the improvement of the perturbative scale errors indeed occurs when higher order 
corrections are included in the calculation of the perturbative terms in the respective 
SCET realization of fl9]). We stress again that the error that we discuss here of 0(Qt/Q) 



9 



due to the approximations in the defining derivation of ([HD in Refs. [SniES] also applies 
to the order- by-order results in Refs. [31]. As we noted above, SCET involves for single 



Z production at the LHC the defining error a/A/M^ ~ 5.7% for the typical hadronic 
transverse size A = 0.3GeV and this is consistent with the approximations made in 
Refs. Such approximations cannot be used for the 1% precision QCD theory that 

we have as our goal for (j3j) here. 

In summary, this last remark is beginning to be manifest as we see in Refs. [27] in the 
comparison between the recent Tevatron and LHC data on the Qt spectra in the 
production and the predictions of ResBos and of Ref. [33]. Indeed, even though a new 
Qr-related variable [36] is used in some of the comparisons, 0* = tan(~(7r — A0)) sin#* = 

p ' T ^ n ^ 4- Q(hr^ where — ( p l _ 2 j s the azimuthal angle between the two 

eptons which have transverse momenta pl T , i = 1,2, and 6* is the scattering angle 
of the dilepton system relative to the beam direction when one boosts to the frame 
along the beam direction such that the leptons are back to back, one sees that these </>*- 
comparisons also show the underlying physical precision error associated with the defining 
formula in as expected, the comparisons are somewhat better than the Qt spectra 
comparisons because this 4>* quantity is more inclusive - two different values of Qt with 
correspondingly compensating differences in the attendant <f>i can have the same value 
of </>*. We want to encourage the LHC experimentalists to continue to produce the real 
Qt spectra in the Z/7* production because we believe, even if there are some irreducible 
experimental systematic errors, these spectra can be very useful in determining which 
theoretical approach is actually correct as we argue in the next section. In particular, 
given the intrinsic physical precision error in fl9]), one has to account for that error when 
the data are normalized to any prediction which is based on the formula therein. We 
would like to stress this latter point. We await the more complete data sets accordingly. 

From the discussions just completed, we see that, in order to have a strict control 
on the theoretical precision in (jTJ, we need both the resummation of the reduced cross 
section and that of the attendant evolution of the {Fj}. We turn now to the latter. 

More specifically, we apply the QCD restriction of the formula in (JJJ) to the calculation 
of the kernels, Pab, i n the DGLAP-CS theory itself and thereby get an improvement of 
the IR limit of these kernels, an IR-improved DGLAP-CS theory [TIE] in which large IR 
effects are resummed for the kernels themselves. The attendant new resummed kernels, 
-P^b are given in Refs. [7H9]. We reproduce the new kernels here for completeness: 

2 



' + 7 (1 - z)T - f q { lq )5{l - z) 



1-z 

Pg(z) = C F F YFS (j q )e^ 1 + il ~ Zy z^, 
P^(z) = 2C G F Y¥S { lG )e¥o { l^l^G + _L_ (1 _ z ya 

+ \(z 1+1G (l -z) + z(l - z) 1+ ^) - f G { lG )5{\ - z)}, 
P%?(z) = F YFS ( 7 G)e^{z 2 (l - zfo + (1 - z )2 z 7g }i (20) 
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where the superscript "exp" indicates that the kernel has been resummed as predicted 
by Eq. (j4]) when it is restricted to QCD alone. Here Cf{Cg) is the quadratic Casimir 
invariant for the quark(gluon) color representation respectively, and the YFS [19] infrared 
factor is given by 

F YFS (a) = e- CEa /T(l + a) 

where T(w) is Euler's gamma function and Ce = 0.57721566... is Euler's constant. The 
definitions of the respective resummation functions ja, 5a, /a, A = q,G are as follows WM 
B 

a s 4C F x l q . a s C F it 2 1 

li = C F — t = -— , S q = -± + ——(-5- - 5), 

7T po Z 7T O Z 

2 2 1 



lq lq + 1 7? + 2 

a s AC G 7 G a s C G vr 2 1 

7g = C g — t=— — , 0q = — I — (— —J, 

7T Po Z 7T O Z 

f ( 7 ) = IV e -^ G + ? 1 

1 ^ 6C g F yfs ( 7 g) 7g(1 + 7g)(2 + 7 g) (1 + 7 g)(2 + 7 g) 1 ' 

j I f22) 

2(3 + 7 g)(4 + 7g) (2 + 7g)(3 + 7g)(4 + 7g)' V ' 

where T(w) is Euler's gamma function and C is Euler's constant. We use a one-loop 
formula for a s (Q) in ( 122]) . so that 

A, = 11 - = 4/3 x , 

where = A^i? is the number of active quark flavors. These new kernels provide us with 
a new resummed scheme for the parton density functions (PDF's) and the reduced cross 
section with the same value of a in (Tj[|): 

F h a ^ F;, a' for 
P G? (*)^ Peseta. 1 ) 

This new scheme has improved MC stability as discussed in Refs. jH] in the attendant 
parton shower MC based on the new kernels there is no need for an IR cut-off '/co' param- 
eter. It is important to note that, while the degrees of freedom below the IR cut-offs in 
the usual showers are dropped in those showers, in the showers in HERWIRI1.031, as one 
can see from @, these degrees of freedom are included in the calculation and are inte- 
grated over in the process of generating the Gribov-Lipatov exponents 7^ in ([2"0]) . We note 
also that the new kernels agree with the usual kernels at 0(a s ) as the differences between 
them start in 0(a 2 s ). This means that the NLO matching formulas in the MC@NLO and 



5 The improvement in Eq. (|20|) should be distinguished from the resummation in parton density evo- 
lution for the "z — > 0" Regge regime - see for example Refs. 35,37 . This latter improvement must also 
be taken into account for precision LHC predictions. 
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POWHEG frameworks apply directly to the new kernels for the realization of exact NLO 
ME/shower matching. 

In the interest of pedagogy, in Fig. 1 we illustrate the basic physical idea, already 
discussed by Bloch and Nordsieck |38j, which underlies the new kernels: an accelerated 



Soft Gluon Cloud 
Gi(£i) * Gt^t) 

JLJL 



q q(l - z) 

q -> q(l - z) + G®Gi ■ ■ ■ <g> G e , I = 0, • • • ,00 



Figure 1: Bloch-Nordsieck soft quanta for an accelerated charge. 

charge generates a coherent state of very soft massless quanta of the respective gauge 
field so that one cannot know which of the infinity of possible states one has made in 
the splitting process q(l) —> q(l — z) + G <g> G\ • • • ® Ge, £ = 0, • • • , 00 illustrated in 
Fig. 1. In the new kernels this effect is taken into account by resumming the terms 

O ((a s ln(^) ln(l - z)) n ^ when z 1 is the IR limit. From ([23]) and ((D), we see that 

when the usual kernels are used these terms are generated order-by-order in the solution 
for the cross section a in ([1]) so that our resumming them enhances the convergence of the 
representation in ([1]) for a given order of exactness in the input perturbative components 
therein. In the next Section, we illustrate this last remark in the context of the comparison 
of recent LHC data to NLO parton shower/matrix element matched predictions. 

3 Interplay of IR-Improved DGLAP-CS Theory and 
NLO Shower/ME Precision: Comparison with LHC 
Data 

In the new MC HERWIRI 1.031 [9] we have the first realization of the new IR-improved 
kernels in the HERWIG6.5 [10] environment. Here, using recent LHC data as our baseline, 
we compare it with HERWIG6.510, both with and without the MC@NLO [llj exact 0{a 8 ) 
correction to illustrate the interplay between the attendant precision in NLO ME matched 
parton shower MC's and the new IR-improvement for the kernels. 

More specifically, in Fig. [2] in panel (a) we show for the single production at the 
LHC the comparison between the CMS rapidity data [39J and the MC theory predictions 
and in panel (b) in the same figure we show the analogous comparison with the ATLAS Pt 
data. Here, the rapidity data are the combined e + e~ — results and the pt data are 
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those for the bare e + e~ case; for, the theoretical framework of our simulations corresponds 
to these data. We do not as yet have complete realization of all the corrections involved 
in the other ATLAS data in Ref. [40J. These results are better appreciated if they are 




Figure 2: Comparison with LHC data: (a), CMS rapidity data on (Z/j*) produc- 
tion to e + e~, pairs, the circular dots are the data, the green(blue) lines are 
HERWIG6.510(HERWIRI1.031); (b), ATLAS p T spectrum data on {Z/'f) produc- 
tion to (bare) e + e~ pairs, the circular dots are the data, the blue(green) lines are 
HERWIRI1.031(HERWIG6.510). In both (a) and (b) the blue(green) squares are 
MC@NLO/HERWIRI1.031(HERWIG6.510(PTRMS = 2.2GeV)). In (b), the green trian- 
gles are MC@NLO/HERWIG6.510(PTRMS =0). These are otherwise untuned theoretical 
results. 

considered from the perspective of our analysis in Ref. [9] of the FNAL data on the single 
Z/j* production in pp collisions at 1.96 TeV. 

More precisely, we direct the reader to the results in Fig. 11 of the second paper in 
Ref. [9]. In that figure, we showed that, when the intrinsic rms pt parameter PTRMS 
is set to in HERWIG6.5, the MC@NLO/HERWIG6.510 simulations give a good fit to 
the CDF rapidity distribution data jl2] therein but they do not give a satisfactory fit 
to the DO pt distribution data [13] therein. In contrast, the corresponding simulations 
for MC@NLO/HERWIRI1.031 give good fits to both sets of data with the PTRMS = 0. 
Here PTRMS corresponds to rms value for an intrinsic Gaussian distribution in px- The 
authors of HERWIG [H] already have observed that, to get good fits to both sets of data, 
one may set PTRMS = 2 GeV. Accordingly, in analyzing the new LHC data, we have set 
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PTRMS = 2.2GeV in our HERWIG6.510 simulations while we continue to set PRTMS=0 
in our HERWIRI simulations. 

We turn now with this perspective to the results in Fig. |2l where we see a confir- 
mation of the finding of the HERWIG authors. One needs to set PTRMS = 2GeV [44] 
in the MC@NLO/HERWIG6510 simulations to get a good fit to both the CMS rapidity 
data and the ATLAS p T data. We again see that the MC@NLO/HERWIRI1.031 sim- 
ulations with PTRMS = at LHC give a good fit to the data for both the rapidity 
and the px spectra. Quantitatively, we use the x 2 /d-o.f. as a measure of the good- 
ness of the respective fits. From the results in Fig. [2] we compute that the x 2 /d.o.f. 
for the rapidity data and the x 2 /d.o.f. for the px data are (.72, .72) ((.70,1.37)) for the 
MC@NLO/HERWIRI1.031(MC@NLO/HERWIG 6510(PTRMS=2.2GeV)) simulations. 
The corresponding results are (.70,2.23) for the MC@NLO/HERWIG6510(PTRMS=0) 
simulations. 

To reproduce the LHC data on the px distribution of the Z/j* in the pp collision the 
usual DGLAP-CS kernels require the introduction of a hard intrinsic Gaussian distribution 
in px inside the proton whereas the IR-improved kernels give in fact a better fit to the data 
without the introduction of such a hard intrinsic component to the motion of the proton's 
constituents. The hardness of this intrinsic PTRMS is the issue as this quality of it is 
entirely ad hoc; it is in disagreement with the results of all successful models of the proton 
wave-function [35], wherein the scale of the corresponding intrinsic PTRMS is found to 
be < 0.4GeV. More significantly, it contradicts the well-known experimental observation 
of precocious Bjorken scaling [45|I4T]; for, the famous SLAC-MIT experiments on the deep 
inelastic electron-proton scattering process show that Bjorken scaling occurs already at 
Q 2 = 1+ GeV 2 for Q 2 = —q 2 with q the 4-momentum transfer from the electron to the 
proton. If the proton constituents really had a Gaussian intrinsic px distribution with 
PTRMS = 2GeV, these pioneering SLAC-MIT observations would not be possible. What 
we advocate now is that the ad hoc "hardness" of the PTRMS = 2.2GeV value is really 
just a phenomenological representation of the more fundamental dynamics described by 
the IR-improved DGLAP-CS theory. This raises the following question: "Is possible to 
tell the difference between the two representations of the data in Fig. |2p" 

One expects physically that more detailed observations should be able to distinguish 
the two representations of the data in Fig. [2J In this connection, in Fig. El we show for the 
Z/Y mass spectrum the MC@NLO/HERWIRIl. 031 (blue squares) and MC@NLO/HER- 
WIG6510(PTRMS=2.2GeV) (green squares) predictions when the decay lepton pairs sat- 
isfy the LHC type requirement that their transverse momenta {p T ,p T } exceed 20 GeV. 
From the results in Fig. [31 wherein the peaks differ by 2.2% for example, we see that 
the high precision data such as the LHC ATLAS and CMS experiments will have (each 
already has over 5 x 10 6 lepton pairs) would allow one to distinguish between the two sets 
of theoretical predictions. 

Continuing in this direction with the discussion, we see that the main differences 
between the three predictions in Fig. [2] (b) occur in the regime below 10 GeV/c. To 
better probe this latter regime, we make a more detailed snap-shot of it in which we 
plot in Fig. 4 the three respective featured theory predictions with the finer binning of 
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Figure 3: Normalized vector boson mass spectrum at the LHC for p T (lepton) > 20 GeV. 



0.5GeV/c instead of the 3.0GeV/c binning used in Fig. [2] (b). The results in Fig. 4 show 
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Figure 4: Normalized vector boson spectrum at the LHC for the ATLAS cuts as 
exhibited in Fig. [2]for the same conventions on the notation for the theoretical results with 
the vector boson p T < 10 GeV to illustrate the differences between the three predictions. 

that the three theoretical predictions have significant differences in the shapes that are 
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testable with the precise data that will be available to the ATLAS and CMS experiments. 
We would again encourage experimentalists to pursue the measurements of both p? and 
0* spectra as these will both be very useful in establishing the correct theoretical approach 
to the respective LHC observations. We will pursue elsewhere (12] other such detailed 
observations that may also reveal the differences between the two descriptions of parton 
shower physics. In this connection especially, we continue to await the release of the entire 
data sets from ATLAS and CMS. 

4 Conclusions 

What we have shown is the following. The realization of IR-improved DGLAP-CS theory 
in HERWIRI1.031, when used in the MC@NLO/HERWIRI1.031 exact 0{a a ) ME matched 
parton shower framework, affords one the opportunity to explain, on an event-by-event 
basis, both the rapidity and the px spectra of the Z/7* in pp collisions in the recent LHC 
data from CMS and ATLAS, respectively, without the need of an unexpectedly hard in- 
trinsic Gaussian p? distribution with rms value of PTRMS = 2 GeV in the proton's wave 
function. Our view is that this can be interpreted as providing a rigorous basis for the 
phenomenological correctness of such unexpectedly hard distributions insofar as describ- 
ing these data using the usual unimproved DGLAP-CS showers is concerned. Accordingly, 
we have proposed that comparison of other distributions such as the invariant mass dis- 
tribution with the appropriate cuts and the more detailed Z/Y px spectra in the regime 
below lO.OGeV be used to differentiate between these phenomenological representations 
of parton shower physics in MC@NLO/HERWIG6510 and the fundamental description 
of the parton shower physics in MC@NLO/HERWIRI1.031. We have further emphasized 
that the precociousness of Bjorken scaling argues against the fundamental correctness 
of the hard scale intrinsic ansatz with the unexpectedly hard value of PTRMS = 2 
GeV, as do the successful models jl5] of the proton's wave function, which would predict 
this value to be < 0.4GeV. As an added bonus, we have pointed-out that the funda- 
mental description in MC@NLO/HERWIRI1.031 can be systematically improved to the 
NNLO parton shower/ME matched level - a level which we anticipate is a key ingredient 
in achieving the (sub-)l% precision tag for such processes as single heavy gauge boson 
production at the LHC. 

Evidently, relative to what one could achieve from the fundamental representation of 
the corresponding physics via IR-improved DGLAP-CS theory as it is realized in HER- 
WIRI1.031 when employed in MC@NLO/HERWIRI1.031 simulations, the use of ad hoc 
hard scales in models would compromise any discussion of the attendant theoretical pre- 
cision. We are pursuing additional cross checks of the MC@NLO/HERWIRI1.031 simu- 
lations against the LHC data. 

In our discussion, we have also spent some amount of time discussing alternative 
approaches to the type of resummation embodied in the IR-improved DGLAP-CS theory; 
for, some of these approaches are in wide use. What we conclude is that the physical 
precisions of these approaches(see Sect. 2), which are based on Refs. are above 
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the 1% precision tag that we aspire, even though there is no contradiction between our 
exact approach and these more approximate methods. 

In closing, two of us (A.M. and B.F.L.W.) thank Prof. Ignatios Antoniadis for the 
support and kind hospitality of the CERN TH Unit while part of this work was completed. 
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